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Abstract 



In Gentile statistics the maximum occupation number can take on unrestricted integers: 
1 < n < oo. It is usually believed that Gentile statistics will reduce to Bose-Einstein statistics 
when n equals the total number of particles in the system N . In this paper, we will show that 
this statement is valid only when the fugacity z < 1; nevertheless, if z > 1 the Bose-Einstein case 
is not recovered from Gentile statistics as n goes to N. Attention is also concentrated on the 
contribution of the ground state which was ignored in related literature. The thermodynamic 
behavior of a ^-dimensional Gentile ideal gas of particle of dispersion E = , where v and 
s are arbitrary, is analyzed in detail. Moreover, we provide an alternative derivation of the 
partition function for Gentile statistics. 



CO 
> 

^ '■ I. Introduction 
o 

— h Quantum statistics is classified as either Bose-Einstein or Fermi-Dirac statistics. As a generalization, 
^ Gentile [1] proposed a kind of intermediate statistics, in which the maximum number of particles in 
"j^ . any quantum state is neither 1 nor oo, but equal to a finite number n. A particularly interesting 

3 ■ case of Gentile statistics is n — > N, where N is the total number of particles in the system. Some 

£3 authors, e.g., Ref. [2], argue that in this case there is no difference between Bose-Einstein and Gentile 
733 statistics. However, this statement is valid only when the fugacity z < 1. In the present paper, we 

g ■ will discuss the case of n — > N with z > 1 and show that in this case the Bose-Einstein case can not 

O ■ be recovered from Gentile statistics. 

>- ! The thermodynamic behavior of an ideal gas obeying Gentile statistics has been discussed in some 
^ | literature [U [3]. Nevertheless, in these works the contribution from the ground state was ignored. 
It is known that such a contribution can be safely ignored in fermion systems; however, in boson 
systems it plays an important role especially in the case of low temperatures and high densities. In 
Gentile statistics, the maximum occupation number n can take on unrestricted integers: 1 < n < oo. 
When n is large, especially n — > N , we will show that the contribution from the ground state 
becomes very important. 

In the present paper, we will calculate the thermodynamic quantities in ^-dimensional space with 
the dispersion relation E = where v and s are arbitrary. A special case of v = 3 and s = 2, 
without the ground state contribution, has been discussed in Ref. [HE]. Our result shows that, the 
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thermodynamic behavior of Gentile statistics, especially the contribution of the ground state for the 
case v > s, depends sensitively on the maximum occupation number n. 

In Sec. II, we first provide an alternative way to calculate the partition function for Gentile 
statistics. Special attention is concentrated on the case of n — > N with z > 1. In this case, Gentile 
statistics will not reduce to Bose-Einstein statistics. In Sec. Ill, we calculate the thermodynamic 
quantities in which the contribution from the ground state is reckoned in. In Sec. IV, we discuss 
the low temperature and high density result for the internal energy, the specific heat, the magnetic 
susceptibility and other thermodynamic quantities. In this case, the ground state contribution plays 
an important role. A summary and conclusions are given in Sec. V. 



II. The partition function and the case of n — >■ N with z > 1 
in Gentile statistics 

For completeness, we first calculate the partition function for Gentile statistics in a way which is 
somewhat different from that in literature. The approach we shall employ is based on multinomial 
theorem. 

Starting with counting the number of states, we express the grand partition function as follows: 

I £ at 

Here 

~ £ = £!^e- (a+feW , (2) 

and Qg is the number of quantum states of ag particles distributed in one level eg (the £-th level) with 
a degeneracy ug. In this occasion, one state can be occupied by at most n particles. It means 

^ = E---E^ (3) 

i=0 

where ctj (i — 0, 1, • • • , n) denotes the number of states which contains i particles, and the summation 
runs over all possible values. The set of {aj} must satisfy the conditions 

OLi = Vi, ia i = a t- ( 4 ) 

i i 

To obtain the grand partition function, we need the multinomial theorem [3]: The coefficient Q 
of x a (0 < a < nu) in the expansion of (1 + x + x 2 + • • • + x n Y is 

i=0 

n n 

where ctj (i — 0, 1, • • • , n) satisfy a i — i a i — a i an d the summation runs over all possible 

i=0 i=0 

values. 

Based on such a relation, we can obtain the analytic form of the H^. The fig in Eq.([3]) (constrained 
by Eq.(jU)) is the coefficient of the term e -( a +feW j n the expansion of 



1 + g-(a+/3££) + e ~2(a+)3e e ) _^ ^ e -n(a+/3e e ) 



(5) 
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As a result, Eq.(T2]) can be calculated exactly: 

~ e -(n+l){a+Pee) _ ]_ 
" l ~ e -(a+(Se e ) _ I 

The result agrees with that found by Gentile []]. 

Obviously, when n = oo the necessary condition on obtaining Eq.Q from Eq.Q is the fugacity 
z < 1 (z = e~ a ); otherwise, will diverge. This means that Gentile statistics will return to Bose 
case only when z < 1. Many authors, e.g., Ref. [2], however, argue that it does not make any 
difference whether one take n = N or n = oo; in other words, Bose-Einstein statistics is recovered 
when n = N. Nevertheless, Eq.(jHJ) will hold so long as n ^ oo even for n ^> 1 and z > 1. The 
behavior of the statistics corresponding to z > 1 and n ^> 1 will be very different from Bose-Einstein 
statistics. In Sec. IV, we will discuss such statistics in detail. 
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III. Thermodynamics with the contribution of the ground 
state 

In this section, we will discuss the thermodynamics corresponding to Gentile statistics in detail. 
Specifically, we will study the contribution of the ground state to the thermodynamic quantities, 
which is ignored in Ref. [TJ [3] . 

The grand partition function for Gentile statistics is 



n 



z n+l e ~(n+l)/3e e 



ze 



-fat 



(7) 



Replacing the summation over I by the corresponding integral and carrying out the integral by parts, 
we have 



lnS 



l 

1 2gV 

x u ^r(|) Jo 



n + l 



£-d£ + lm 



1 - z 



n+l 



(8) 



while we have used the asymptotic formula in v dimensions with the dispersion relation E 

V r2W 2 l 



2m" 



Y h u J r(|) s v ' 



and introduced e = kT£. It should be emphasized that, the parameter m in the dispersion relation 
can be explained as the mass of particles only when s — 2. Here A is the thermal wavelength 



A 



h s 



l/s 



27r s / 2 mkT 



(9) 



and g is a weight factor that arises from the internal structure of the particles (the number of internal 
degrees of freedom). We have split off the term in Eq.® corresponding to e = 0, which we explicitly 
retain because, when v > s, the density of states is zero. The physical interpretation of such a term 
is similar to the Bose case: Sometimes the single term e = may be as important as the entire sum. 
The term e = represents the contribution of the ground state of the system. The contribution of a 
state is proportional to the number of particles which are accommodated in the state. As we know, 
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particles tend to occupy the lower state as possible as they can. At the low temperature limit the 
particles all want to occupy the ground state; thus the particle number of the ground state nearly 
equals to the maximum occupation number of the state. Therefore, when T — > the contribution 
of the ground state will be proportional to the maximum occupation number n. For fermions, the 
contribution from the ground state can be ignored since n = 1; however, for bosons, n = oo, the 
contribution from the ground state becomes very important in the low-temperature limit. Indeed, 

when v > s, the term e = diverges in the Bose case and causes the phase transition Bose-Einstein 

condensation. In Gentile statistics the maximum occupation number n is finite, so the term e = 
does not diverge and there are no phase transitions like Bose-Einstein condensation. However, when 
n is large the influence of such a term may be very important. In the following we will show this 
influence on the thermodynamic quantities. Moreover, to obtain the result corresponding to v < s, 
we only need to drop the term e = 0. 

In the theory of Bose-Einstein and Fermi-Dirac systems we come across integrals of the type 

I rod ^ 1 f°° X a ^ 

9o{z) = t^-t / -zr~ x T dx and fM = FTT / -l x i 1 dx - 

1 [a) Jo z 1 e x — 1 1 (a) Jo z 1 e x + 1 

To deal with the systems with the grand potential Eq.(jSJ), we can introduce the function 



K{z) 



Via) 



n+1 



x 



~ x dx. (10) 



z — lgx -y ~— (n+l)g(n+l)ic \ 

The function h a (z) will return to Bose-Einstein or Fermi-Dirac integrals when n equals to 1 or 
oo, respectively. A simple process of differentiation brings out the relationship between h a (z) and 
ha-x(z): 

d 

h a -i(z) = z—h a (z). (11) 
oz 

In terms of the function h a (z) we can express the grand potential Eq.flHJ) in the following form: 

1 laV v z n+1 — 1 

1- = ElnS- = T ^V ( - + l )ki+lW + .„ — . (12) 

The last term describes the contribution from the ground state. For bosons, this term returns to 
— ln(l — z) with n — > oo and z < 1, which is in connection with the Bose-Einstein condensation; for 
fermions, this term may be ignored since the ground state of a fermion system can contain only one 
particle. In Gentile statistics 1 < n < oo, the contribution of the ground state depends on the value 
of n. When n is close to 1, like the Fermi case, the contribution of the ground state can be ignored; 
however, for v > s, when n is comparable in magnitude to the total number of particles, the ground 
state contribution becomes important. Notice that, although there are no divergences, consequently 
there are no phase transitions, the influence which is related to n also can be very large. Ideally, 
phase transitions display only when infinities appear. In practice, however, infinity only means very 
large, so when n —>■ N, though there are no divergences, the behavior of the system will be very 
similar to phase transitions. 

Since the grand potential Eq. ffT2l involves the contribution corresponding to e = 0, from here 
onward, the thermodynamic quantities which reckon in the influence of the ground state follow 
straightforwardly. 

Eliminating z between the two equations 

P 1 2 ^ + 1) ^ 1 z n+1 - 1 

kf = v 111 a = ^^fW h ^ {z) + v ln ~^T> 
v = v = v z d-z^ = ^^rtlT^ ^ 
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we obtain the equation of state of ideal Gentile gases. Here, the last terms of these two equations 
describe the contribution from the ground state, and 



N = n - (- 



1 



71+1 



(14) 



z - 1 z n+1 - 1 

is the average occupation number of the ground state. Easily to see that when n = oo and z < 1 
the occupation number iVo = jf^. This is just the case of bosons. In the case of fermions, Nq = 
and its contribution is negligible compared to the first term in any cases. In Gentile statistics 
1 < n < oo, the influence of such a term depends on the value of n. When n is large this influence 
becomes important. 

The internal energy of the ideal gas is given by 



— -— — " 



N °) u kT k ^ Z) 
Ns ha(z) 



For the Helmholtz free energy of the gas, we get 

F = Nfi- PV = NkT 
and for the entropy 



T 



1 N } hvJz) 

s v ' 

hE +1 (z) 



n+l 



N z 



- 1 



N s hz{z) 



1 , z n+1 - 1 

In z H In 

N z-l 



The specific heat capacity is given by 



(15) 



(16) 



(17) 



dU 
df 



V -Nk(l - ^ 

s N ' 



1 + 



,v_ ^+i(f) 
j %(z) 



iV-iVo V(^-i) 2 



In the derivation of CV we have used the relation 



[n + 1) V +1 




1 Mf) 
/iii-i(^) 



^2; 



1 + 



1 



(n + l) 



n -1 



N — Nn 



I) 2 



(19) 



2 / h^{z) 

f / s \ / . 

The existence of ground state term 



which can be obtained from the expression of N in Eq. (Tl3|) 
makes the formulas somewhat complex. 

Comparing with the thermodynamic quantities given by Ref. [TJ [3], we can see that there is an 
extra factor 1 — ^ in our result. The addition of such a factor to the thermodynamic quantities 
brings in important information: the influence of the ground state. If we ignore the contribution 
from the ground state or, equivalently, take iVo = 0, the above thermodynamic quantities will return 
to the results given by Ref. [H [3]. However, since in Gentile statistics N Q can take any value, the 
influence of the ground state must be reckoned in. For instance, when z ^> 1 and iVo w iV so that 
the factor 1 — ^ — > 0, the thermodynamic quantities, such as the internal energy, will approach 0. 
The ratio jt- reflects the difference of the various intermediate statistics corresponding to different 
n. Notice that, the factor 1 — ^ comes from the term corresponding to e = 0. Such a term is needed 
only when v > s. 

For studying the properties of thermodynamic quantities in detail, we first need to discuss the 
behavior of z as determined by the second equation of Eq. (fl3l) . namely: 
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Instead of the Bose-Einstein and the Fermi-Dirac integrals g»{z) and fz(z), there is a function ha.(z). 
To determine the range of the fugacity z, we rewrite the factor in function ha.{z) in the following 
form: 



£ - 1 f "+1 - 1 



1 n + 1 



f n-i + 2£ n " 2 + h n 

£n _|_ _| 1_ i 



(21) 



where £ = z~ 1 e a: . It is obvious that for any real values of z, hz.(z) is a bounded, positive function of 
z. 

From Eq. fl2ip . one can learn that if we expect ha(z) to return to the Bose-Einstein integral the 
necessary condition z < 1 is needed so that ha(z) is bounded. In other words, if we expect the 
above formulas to return to the Bose case when n = oo, we need an additional restriction z < 1. For 
comparison we recall that < z < oo in Fermi-Dirac statistics and Gentile statistics. 

IV. Low temperatures and high densities 

In this section, we will focus on the behavior of the statistics corresponding to 2; > 1. For \ v jv ^> 1 
the average thermal wavelength is much greater than the average interparticle separation. It is easy 
to see that this is equivalent to the requirement z ^> 1. It can be expected that in this case the 
contribution from the ground state will become important. It is noteworthy that, in this case, the 
result will not return to the Bose case even when n — > N. To obtain Bose-Einstein statistics from 
Gentile statistics, one needs to not only perform the limit n — > N, but also restrict the fugacity 
z < 1. Thus the asymptotic results of the thermodynamic functions for low temperatures and high 
densities given in this section will not return to the Bose case since the asymptotic expansion is based 
on z > 1. (We need not to consider the case of n — > iV with z < 1 because this is just the Bose case 



To study the behavior of h a (z), which contains two integrals, for large z, we introduce the variable 
t = lnz so that the first integral in Eq. fllQI) can be expanded in powers of t: 



while the lower limit t in one of the integral is approximately replaced by 00. By the same procedure, 
we can derive the second integral: 



0) 




(22) 




ti +L / u \ 

( - + 1) FT +a ^0) 



1 



r(i + j)c(i+i). 



(23) 



(n + 1)3 
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Then we obtain: 



h„(z) 



+ 45n (1 



-t a n 



1 + 



1 



7T 



1 



3 n + 1 



-en a 



i)r 



(n + 1) 3 

Substituting Eq.(f24]) into Eq.(JT3D, we obtain 



)a{a - l)(tr - 2)(<7 - 3)t -4 + 



(24) 



iV 

r 



29 



45 n 



-t*n 



1 + 



7T 



3 n + 1 s s 



(--i)r 



,v 



V 



(n + 1 r s s 



-(--l)(--2)(--3)t" 4 + 



(25) 



To obtain the thermodynamic functions for low temperatures and high densities we first write 
down the expansion for the chemical potential from Eq. fl2"5l) : 



u = e F 



The expansion parameter is — . This expression implies that when v = s the chemical potential 



7T 2 1 ,v . ,kT . 

- 1) — + ■ 

3 n + Vs )K e F J 



(26) 



/i is equal to the Fermi energy e F at an arbitrary temperature. A special case is free particles in 
two dimensions, in which v = s = 2. The relation between e F and the Fermi energy e F 



fermion 



111 



Fermi-Dirac statistics is 



e F — e F 



fermion 



/N — n\~ 

V Nn 



where 



^fermion 



11 

2m 



2^ u N 
g y 2 V 



(27) 



(2f 



In other words, e F can be regarded as an analogue of the Fermi energy in Gentile statistics, and e F 

• ii , i fermion i -i 

will return to e F when n — 1. 

With the help of Eq. ([2"6"j) . the asymptotic expansion of the internal energy is then given by 

U N-n 1 



N 

The pressure of the gas is: 



N s. + i 

V 



1 + T ^t(- + 1)(— ) 2 + 
3 n + 1 s e F 



(29) 



1 N-n 

fTT~ 



7T 



1 + 
1 



3 n + Vs ){ e F ; 



) 2 + 



3 n + 1 s 

The free energy of the system follows directly from Eq. (1291) and Eq. (1501) 
F PV N-n 1 



iv =/i 



N 



N i + 1 



1 (- + 1)( ) + 



(30) 



(31) 



whence we obtain the entropy of the system: 



S N-n27r 2 1 vkT 
~Nk ~ ~lv 3~n+i77^ + '"' ^ 
From Eq. fl29l) . the specific heat of the gas can be obtained 

C v _ N - n 2tt 2 1 v kT 

M~^V 3~n + ls e]T + '"" ^ 
It is easy to see that there is a factor (N — n) appearing in each expression of the thermodynamic 
quantities, where N is the total number of particles and n is the maximum occupation number of a 
state in Gentile statistics. When n is of the order of N, the values of the thermodynamic quantities 
will be strongly suppressed by this factor. It is worth comparing Cy to the specific heat of the 
ideal Fermi gas (jff rmion so as to discuss the difference between Gentile statistics and Fermi-Dirac 
statistics: 

C v = Vl C v ermion , (34) 

where 

"■ = (—) ~i (35) 

Similarly, the magnetic susceptibility per unit volume of a system is given by 

X = V2X fermton , (36) 



where 



and y^ erm%on is the magnetic susceptibility of a fermion system: 

„ 0^/2 

fermion _ V * n (n^S^ (Jermum\Z-l 



X } ~ = ^f^(2m)t( e ^~)7-\ ( 3J 



Obviously, the influence of the ground state depends sensitively on the maximum occupation number 
n when v > s. 



V. Conclusions and outlook 

In summary, we first discuss the case n — > N. In some literature, e.g., Ref. [2], the authors argue 
that such a case is just the Bose case. In this paper we point out that this argument is valid only 
when the fugacity z < 1. When z > 1, Gentile statistics with n — > iV does not return to Bose- 
Einstein statistics and in such a case the contribution from ground state plays a very important 
role. Moreover, we discuss the contribution of the ground state in Gentile statistics which is ignored 
in the previous literature. The result shows that the ground state contribution becomes important 
when the maximum occupation number n is large, especially in the case of low temperatures and 
high densities. The results given by current literature focus only on the 3-dimensional space with 

2 

the dispersion relation E — In the present paper, we give the thermodynamic quantities for 
an arbitrary dispersion relation E = ^ in arbitrary v- dimensional space and find that the result 
depends sensitively on the maximum occupation number n when v > s. 
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As a generalization of Bose-Einstein and Fermi-Dirac statistics, fractional statistics has been 
discussed for many years [5j [6] . The particles obeying fractional statistics must not be real particles. 
As an effective method, however, the introducing of such imaginary particles can be used to deal 
with some complex interaction systems [7J. Generally speaking, there are two ways to achieve the 
fractional statistics: (1) One is based on counting the number of many-body quantum states, i.e., 
generalizing the Pauli exclusion principle [H El E] • The most direct generalization is to allow more 
than one particles occupying one state, and such an approach gives Gentile statistics [Ij. Another 
way is to introduce a parameter a valued from to 1 in the expression of the number of quantum 
states Bose-Einstein or Fermi-Dirac statistics becomes its limit case when a is or 1, respectively. 
(2) Alternatively, the fractional statistics can be achieved by analyzing the symmetry properties of 
the wave function. As is well known, the wave function will change a phase factor when two identical 
particles exchanges. The phase factor can be +1 or —1 related to bosons or fermions, respectively. 
When this result is generalized to an arbitrary phase factor e 10 , the concept of anyon is obtained 
[6]. Comparing these two approaches to fractional statistics a question arises naturally: Are there 
any connections between them? To answer this question we need to build a bridge between the 
generalized Pauli principle (in which the maximum occupation number is extended to an arbitrary 
integer n) and the exchange symmetry of identical particles. We will discuss this problem elsewhere 

m. 
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